Intrinsic geometry approach to surface kinetic roughening 
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A model for kinetic roughening of one-dimensional interfaces is presented within an intrinsic 
geometry framework that is free from the standard small-slope and no-overhang approximations. 
The model is meant to probe the consequences of the latter on the Kardar-Parisi-Zhang (KPZ) 
description of non-conserved, irreversible growth. Thus, growth always occurs along the local normal 
direction to the interface, with a rate that is subject to fluctuations and depends on the local 
curvature. Adaptive numerical techniques have been designed that are specially suited to the study 
of fractal morphologies and can support interfaces with large slopes and overhangs. Interface self- 
intersections are detected, and the ensuing cavities removed. After appropriate generalization of 
observables such as the global and local surface roughness functions, the interface scaling is seen in 
our simulations to be of the Family- Vicsek type for arbitrary curvature dependence of the growth 
rate, KPZ scaling appearing for large sytems sizes and sufficiently large noise amplitudes. 

PACS numbers: 68.35.Ct, 05.10.Gg, 64.60.Ht, 81.15.-z 



I. INTRODUCTION 

Kinetic roughening in non-equilibrium intcrfacial 
growth has been investigated for a long time both as 
a fundamental problem in statistical mechanics and as 
a model to relevant physical systems, such as dynamics 
of thin films, fluid flow, flame front propagation, or un- 
differentiated biological growth Starting from a flat 
shape, for many interfaces it is generally observed that 
they roughen until eventually a steady state is attained, 
in which the interface roughness, W (root-mean-squarc 
deviation around the mean interface position) scales in 
a nontrivial way with the system size L, as W ~ L a , 
where a is the so-called roughness exponent. At short 
times in these systems, the roughness scales as a power 
law of time, W(t) ~ t 13 , with /3 being called the growth 
exponent [l[ . The saturation roughness can be measured 
at smaller length scales I than the system size, another 
power law ensuing, W(l) ~ l aioc . If the local and global 
roughness exponents are equal, ai oc = a, the system fol- 
lows the so-called Family- Vicsek (FV) Ansatz, which is 
a generalization of that expected for equilibrium critical 
dynamics [l[ ; otherwise, anomalous scaling occurs 0, Q 
in which interface fluctuations follow different laws at 
small and large scales, a fact that is not uncommon when 
the interface profile displays many large jumps and high 
slope values 0, [f|. In general, the values of the critical 
exponents a, cq oc , and (3 identify universality classes that 
play an important role in the classification of nonequilib- 
rium processes in spatially extended systems Q. 

Fluctuations arise in the above physical systems due 
to a variety of sources, such as non-homogeneous driv- 
ing, randomness of growth events, etc. Dynamically, and 
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in the absence of specific conservation laws that constrain 
the evolution of the interface [2| , they compete with ir- 
reversible growth, which can be simply represented by 
a constant growth rate along the local normal direction 
to the interface, and are attenuated through various re- 
laxation processes. The latter can be differential absorp- 
tion due to curvature effects (e.g., condensation of atoms 
from a vapor phase is more effective at surface valleys 
of a solid than at surface peaks) Qito surface diffusion 
(atoms move from peaks to valleys) [cM10l|. etc. Different 
relaxation mechanisms can give rise to different univer- 
sality classes, some of which may in general appear as 
transients during the system evolution. 

One of the most fruitful approaches to surface kinetic 
roughening is the study of stochastic partial differential 
equations (SPDE) of a certain height function, h(x,t), 
that represents the surface height position at time t above 
point x on a reference plane, in the so-called Monge gauge 
[111 Il2j . This implies the assumption that the interface 
has no overhangs. Moreover, most height equations in 
this context rely on a small-slope approximation that 
allows to strongly simplify the analytical formulation. 
However, this simplification may be at the cost of in- 
accuracies in the physical description. For instance, sur- 
face diffusion is thus taken along the horizontal substrate 
plane (instead of along the interface) and is severely over- 
estimated in regions of high slopes, giving rise to large 
grooves as in the so-called linear Molecular Beam Epitaxy 
(MBE) equation [H[. 

The very same Kardar-Parisi-Zhang (KPZ) equation 
that is one of the most important SPDEs within 
the study of kinetic roughening — and within non- 
equilibrium statistical mechanics indeed — was actually 
formulated within both the no-overhang and the small 
slope approximations. The fact that the corresponding 
roughness exponent «kpz < 1 Q| makes this approxi- 
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mation self-consistent, in the sense that slopes remain 
small on average. Nevertheless, the derivation of this 
equation being based on rather general considerations, 
a natural question arises as to what is the effect of re- 
laxing these approximations. Actually, early reports [l4| 
already claimed that relaxing the no-overhang approxi- 
mation in a KPZ-type equation leads to a larger value for 
the growth exponent than that characteristic of the stan- 
dard KPZ equation, although no information was pro- 
vided on the value of the roughness exponent, nor are 
later confirmations of this result available. One of the 
motivations of the present work is to address a more sys- 
tematic study along these lines. The question becomes 
all the more interesting since experimental realizations of 
KPZ scaling are remarkably scarce, in spite of the wide 
asymptotic universality that it is expected to have on 
theoretical grounds [2, Nevertheless, this univer- 

sality class has been indeed verified in experiments very 
recently [13, including detailed predictions [H, on 
the asymptotic height distribution for the case of growth 
in a circular geometry. Analogous predictions are also 
available for the simpler case of a band geometry [2C} in 
which boundary conditions are periodic or free, and the 
initial interface is a straight, horizontal line. Still, these 
remain to be experimentally assessed [ItJ- 

Actually, the relation between the scaling properties 
of a given SPDE in band and in circular geometries has 
been the focus of a recent controversy that has arisen 
from conflicting analytical predictions |2l| - l23| . These un- 
avoidably require approximations, since even equations 
that are linear in a band geometry become nonlinear in 
the circular setting. Thus, in order to address the issue, 
it would be interesting to be able to study numerically 
the same SPDE in both cases. The circular geometry re- 
quires naturally a multiple valued interface (as seen from 
a Mongc gauge) and arbitrary interface orientations, the 
additional existence of overhangs being again the most 
complete situation to be considered. This constitutes a 
second motivation for the present work, and makes the 
formulation of such an equation in a band geometry — as 
in the seminal formulation by KPZ — a natural prelim- 
inary step. Note, moreover, that the observables (cor- 
relation functions, etc.) also need to be adapted to the 
type of geometry in which the continuum model is for- 
mulated, see e.g. [2~H and references therein. Again from 
this point of view, starting out with the band geometry 
seems a natural choice. 

Previous attempts have been made in order to study 
continuum models of kinetic roughening beyond the two 
approximations mentioned. For instance, Sicgcrt and 
Plischke [25| . Maritan et al. 14 1 or Family and Amar 
[lH considered various relaxation operators still in the 
Monge gauge, but beyond the small slopes approxima- 
tion, thus studying very complex and highly nonlinear 
SPDEs that are hard to deal with even numerically. Ac- 
tually, later developments have questioned the reliability 
of the schemes employed, since the standard finite dif- 
ference discretization of nonlinear terms has been shown 



to lead to numerical inaccuracies and instabilities (see an 
overview in e.g. [26j]). Approximations in the correspond- 
ing SPDEs are also usually performed in the study of con- 
tinuum growth equations in circular geometry, typically 
including the neglect of overhangs as well as additional 
simplifications, see e.g. 0, H3| and references therein. 
The specific case of the KPZ equation has been addressed 
in some of these works [3, [H| and elsewhere, see e.g. 
[28l [29| . although overhangs seem to have been consis- 
tently neglected in all these SPDE-based approaches. 

Also, phase-field studies have been carried out since 
the pioneering work of Kcblinski et al. [3(| HH- How- 
ever, while in this class of studies arbitrary interface 
geometries are available, computational expense is suf- 
ficiently high that, to our knowledge, no complete sys- 
tematic study of scaling properties is available so far. 
Moreover, connection with the KPZ equation is at best 
indirect via the asymptotic behavior, if appropriate. 

In this work, our approach is to study a one- 
dimensional interface as a 2D curve in a band geome- 
try as described above, that evolves under a continuum 
equation whose terms are determined by its intrinsic ge- 
ometry. Actually, these are the same as those that led 
to the formulation of the KPZ equation, but now free 
from the no-overhang and the small-slope approxima- 
tions. We study the ensuing model numerically, carrying 
out simulations with an adaptive number of points, and 
doing measurements in a geometrically natural way. The 
equations do not distinguish between the horizontal (sub- 
strate) and vertical (growth) directions, their difference 
residing only in the initial and boundary conditions. 

The paper is organized as follows. ScctionlTllintroduces 
the physical model we study and its geometric interpre- 
tation. Next, in Section IIII1 we discuss the numerical 
implementation of the simulations. In Section IIVI the 
observables to be measured are considered, while results 
are presented in Section [V] Further discussion of these 
is presented in the final Section PVTI together with some 
conclusions and indications of potential new directions 
for future work. 



II. INTRINSIC GEOMETRY FOR SURFACE 
GROWTH IN ID 

A. Intrinsic Geometric Preliminaries 

Our ID interface is given mathematically by a curve 
7 in a 2D band geometry, that we take to be an infinite 
cylinder of circumference L, i.e., r : [0, 1] S 1 x K. 
Of course, r(0) = r(l). The curve is assumed to wrap 
around the cylinder exactly once, thus separating a solid 
phase (below) from a dilute phase (above). 

At any point on the curve r, let us consider the uni- 
tary tangent vector field, u t , and the normal vector field 
pointing towards the dilute phase, u„. The curvature 
K(r) along the interface will be given by the inverse of 
the radius of the best fitting (osculating) circle. 
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For any scalar field defined on the curve 4>(r), we 
define its 7-gradient as the derivative along the curve, 
parametrized by its arc- length s, pointing along the tan- 
gent: 

V 7 0(r( S )) = ^u t (r( S )). (1) 

The divergence along the curve of such vector field allows 
us to define a Laplacc-Bcltrami operator as 

V?«r<»)) = P> 

B. Physical Model 

In order to have a physical image in mind, let us con- 
sider a system in which atoms are interchanged between 
a dilute (vapor or liquid) and a solid (aggregated) phase, 
in such a way that the latter grows at the expense of the 
former, as in the production of thin films [la ]. We will 
consider absorption and desorption to be non-directional 
phenomena, and the possibility of atomic diffusion at the 
interface. The chemical potential difference at an inter- 
face point, fJ.(r), will be the only driving force for growth. 
On the solid side of the interface (subindcx s), the chem- 
ical potential will be given by Q 

/i s (r) = n S:0 - (j, s>1 K(r), (3) 

where ^ St o and /i Sj i > are empirical parameters, and 
the local curvature K(r) provides a measure of the num- 
ber of local "dangling bonds" . We will consider K to be 
positive when the interface is concave as seen from the 
dilute phase. On the other hand, the chemical potential 
in the dilute ("vapor", subindex v) phase is /z„(r) = (i v .o- 
Given that the solid phase is growing at the expense of 
the dilute phase, we can take ^ Si o < and /j, v .q = with- 
out loss of generality On the other hand, particles 
may move along the surface in their attempt to minimize 
their chemical potential. This process is called surface 
diffusion (32I |33| . Locally, the particle current will be 
proportional to the surface gradient of fi s , 

j s (r) oc -V 7/ u s (r). (4) 

Conservation of the number of particles requires that 
they aggregate to the solid phase at a rate proportional 
to the divergence of this current , thus 

V 7 • j s (r) cx -V 7 /i s (r) = /i M V^A'(r). (5) 

In line with the seminal approach by KPZ Q, a point 
at the interface r is assumed to move along the local 
normal direction, with a velocity that is proportional to 
the rate at which particles accumulate there, namely, 

d t r = (-r(/x s (r) - fi v ) - V 7 ■ j s (r) + r] v (r,t)) u„, (6) 

where T > is a mobility constant Q and r\ v is a zero- 
average noise term — describing, e.g., fluctuations in the 



rate of growth events — that we will take as Gaussian, 
with zero average, and white both in space and time, 

(^(r, t) Vv (r' } t')} oc <5(r - r')S(t - t'). (7) 

Note that the delta function in ([7]) is evaluated at points 
on the interface. Thus, it does not correspond to a delta 
function on the horizontal substrate plane. 

The general form of the ensuing evolution equation is 
thus 

8 t r = (A + A!K(r) - A 2 Vp<(r) + A n?7 (r, t)) u„. (8) 
This equation contains four terms: 

• A constant term, Aq, measuring the average differ- 
ence between the chemical potentials in the dilute 
and solid phases. 

• A term proportional to the curvature, with coeffi- 
cient A\ > 0, having the physical interpretation of 
surface tension Q. 

• A term proportional to the Laplacian of the curva- 
ture field along the curve, with coefficient A2 > 0. 

• A white noise term, due to the fluctuations in 
growth events, with variance A\. 

In this work we study the scaling properties of Eq. (JSj) 
for Ai = 0, leaving the general case for further work. 
Therefore, the equation under study will be 

d t r = (A + iliJf(r) + A n r](r,t)) u„. (9) 

This SPDE is analogous to the KPZ equation in the sense 
that it incorporates the same basic ingredients, namely, 
growth along the local normal direction, relaxation by 
surface tension, and fluctuations, but there are several 
important differences: (a) growth is always normal to the 
interface, while in the KPZ equation this condition holds 
only approximately; (b) we employ the full curvature, 
which is nonlinear if written in the Mongc gauge, and (c) 
noise is normal to the interface, which means that, as seen 
in the Monge gauge, it is multiplicative and correlated. 

In order to represent an interface, the curve should be 
simple for all times, i.e.: no multiple points are allowed. 
Therefore, Eq. <j9j) should be supplemented with a con- 
dition regarding the treatment of self-intersections. Our 
criterion is simply removal, keeping the component which 
wraps around the cylinder. 

Finally, an initial condition is required. Unless oth- 
erwise specified, the interface will start as a horizontal 
line. 



III. NUMERICAL SIMULATIONS 

The numerical discretization of Eq. © has been car- 
ried out keeping in mind the desideratum of geometric 
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naturality, using concepts borrowed from discrete dif- 
ferential form theory [34(. The idea of our numerical 
procedure is to employ an adaptive scheme that docs 
not particularly privilege any system direction, and that 
avoids e.g. standard explicit finite differences discretiza- 
tion of the various differential operators appearing, since 
this requires an explicit parametrization of the interface 
and is prone in this context to numerical instabilities, 
see e.g. [HI and references therein. The guiding physical 
principle is to use discrete counterparts of the operators 
appearing in Eq. (J9J), which are geometrically motivated. 

The interface is simulated numerically by a chain 
of points {Pi, . . . , Pn}, with adaptive N and periodic 
boundary conditions, Pn+i = Pi- It is implemented in 
the computer in the form of a linked list, ensuring easy 
insertion and deletion of points. This set fulfills the fol- 
lowing continuity condition |36f: the distance between 
any couple of neighboring points should always be in the 
interval \l m in > l-max] ■ Let us use the symbol dij for the 
distance between any two points of the set. Then, at 
every time step, this condition is ensured using the fol- 
lowing algorithm for every couple of nearest neighbors 

Pi, Pi+l- 

• If < Imin, remove Pj+i from the list. 



If d 



i,i+l 



> It 



add a new point Pj+i = (P, 



Pj_l_i)/2 (and shift the indices appropriately). 

Application of this algorithm ensures that the represen- 
tation of the interface is homogeneous. It is important, 
in all simulations, to generate the initial points randomly, 
in order to avoid lattice artifacts. 

Next we provide a discrete approximation to the dif- 
ferent geometric fields needed to simulate Eq. ([9]). The 
tangent vector field, u t {x) is approximated by {ut)i,i-\-i ~ 
(P;+i — Pi)/di t i+\. Notice that these values are not at- 
tached to single points, but to links between two con- 
secutive points. With it we can approximate the normal 
vector field, {M n )i.i+i, finding the orthogonal vector to 
(ut)t,j+i, pointing towards the dilute phase. 

The curvature field, K{r), is approximated using the 
triangle Pi_i, Pi, Pj+i- We define K{Pi) as the inverse of 
the radius of the circumscribed circle passing through all 
three points. Let on be the angle at P{. Then, 



Ki = K(Pi) 



2sin(ai) 
di-i. 



(10) 



Notice that, as defined, the sign of the curvature depends 
on the sign of on. As in the previous Section, it will be 
positive when the interface is concave, as seen from the 
dilute phase. 

As our last step, we will determine the approximation 
to the gradient and Laplacian along the curve of any 
scalar field </>(r), using the notation </>,; = 4>{Pi), 



(V, 



H,i+1 



H,i+1 



■(ut) 



(11) 



These values are attached to links, not to points. The 
approximation to the Laplace-Beltrami operator is com- 
puted as the (numerical) divergence of that expression: 



H,i+1 



H,i+1 



(12) 

This formula may be formally justified using the notion 
of Hodge dual [34[ . 

As pointed out before, the interface may tend to create 
self-intersections, leading to cavity formation. Our code 
detects such self-intersections and removes one compo- 
nent, keeping the one which wraps around the cylinder. 
The detection algorithm is one of the most delicate and 
costly ingredients of our numerical scheme. 

Finally, concerning the time evolution, the SPDE (J9j> 
is simulated using the usual Euler-Maruyama algorithm. 
We have ensured that our At is small enough for the com- 
putations by repeating a selection of them with smaller 
values, always obtaining convergence. Regarding the 
noise amplitude, it is important to remark that, in the 
discretized SPDE, it must be multiplied by {As)~ 1 / 2 , 
where As is the arc-length interval corresponding to each 
point. This guarantees that the noise amplitude does not 
depend on the local interface orientation [l2| ■ 

Additional checks have been performed on our algo- 
rithm in order to reproduce expected features in the be- 
havior of Eq. © (not shown): (1) In the deterministic 
case {A n = 0) with Aq = 0, we have performed a linear 
stability analysis of Eq. ^ , plotting the relaxation time 
of a sinusoidal initial condition with wave-vector k, ob- 
taining the expected diffusion-like relation between typ- 
ical relaxation time and spatial frequency, r <~ fc~ 2 . In 
this sense, the deterministic case of Eq. © reproduces 
(non-linear) diffusion. (2) For small non-zero Aq, with 
A n = A\ = 0, a sinusoidal perturbation evolves towards 
the typical deterministic KPZ arcade. 



IV. MEASUREMENTS 

The standard basic observable in kinetic roughening 
theory is the global interface width or surface roughness. 
It is usually defined within the Monge gauge as 



W M 



ongc 



= ({h-hy 



1/2 



(13) 



where h is the height of the interface, as measured from 
a horizontal reference line, bar stands for space average 
across the full system and brackets denote noise real- 
izations. This definition, therefore, privileges the hor- 
izontal orientation. Moreover, if the interface presents 
overhangs, the very use of h{x) becomes ambiguous. 

In many applications and in experimental studies of 
e.g. thin film dynamics [37| . it is customary to remove 
the plane, i.e., finding the optimal orientation which can 
be taken as the horizontal. Our definition of interface 
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width goes beyond this, being completely intrinsic at ev- 
ery length scale. Given any portion of interface, proceed 
in the following way: 

• Find the straight line which makes the best fit to 
the curve. 

• Compute the root mean square distance of the in- 
terface points to this line. 

Skipping the first step and considering the best fitting 
line to be horizontal amounts to using the Mongc gauge 
definition of interface width. Note that the fit is not an 
usual least-squares fit. In such a case, the minimized 
distances are taken along vertical lines. In our case, we 
minimize the actual distances, i.e., along a diagonal line, 
from the curve points to the fitting line. 

Likewise, we define local width w(l,t) at any length 
scale I in the usual way: sample the surface using win- 
dows of size I, and measure the average width obtained, 
following the previous procedure. Of course, it is ex- 
pected that w should be an increasing function of I. At 
the maximum length-scale, / = L, we retrieve the global 
width W(t) = w(L,t). We should remark that the win- 
dows should be random. Otherwise, lattice artifacts are 
bound to appear. 

The Monge-gauge roughness at small scales, 
linij^o w Monge(^ 0: describes the slopes statistics. 
In our case, if the surface is smooth below a certain 
cutoff length, such limit will be close to zero. If the 
interface is defined, as in our numerical simulations, by a 
set of linked line segments whose lengths lie in a certain 
interval, there is an upper bound for the roughness at 
each scale. 

A power-law behavior is expected, w ~ P loc , for length 
scales under a certain correlation length, I < Length 
scales above this value are still uncorrelatcd. As time 
evolves, increases until it reaches the system size. 
At that time all length scales are correlated, the system 
saturates to a stationary state and its width does not 
increase any further. The correlation length grows with 
another power law, £(t) ~ t x l z , defining the dynamic 
exponent, z. Using this definition, it is easy to check 
that the saturation time scales as t s ~ L z [l|. Scaling is 
regular, or FV type, when the roughness of the correlated 
scales stays fixed for all times. In that case, a\ oc = a and 
it can be shown that z = a/ (3 On the other hand, 
in the presence of anomalous scaling, a\ oc ^ a and the 
roughness of the correlated scales keeps increasing after 
correlation, only saturating when = L [1, Q . 

V. RESULTS 

Our main result is that, in spite of indeed featuring 
overhangs and relatively large local orientation angles, 
interfaces described by Eq. © follow a standard Family- 
Vicsek scaling, with exponents that fit into the standard 
universality classes. For smaller systems or low values 



of A n , we obtain exponent values that are not far from 
those of the Edwards- Wilkinson (EW) universality class 
in ID, i.e., a^w = ct EWJoc = 1/2, f3 EW = 1/4, z E w = 2. 
This is the behavior associated with the weak coupling 
regime of the KPZ equation. Indeed, recall that in the 
Monge gauge the (naive) KPZ coupling constant scales 
as the square of the noise amplitude [l[ . If L and A„ are 
large enough, and within error bars, we obtain exponent 
values compatible with those of the (strong coupling) ID 
Kardar-Parisi-Zhang universality class, namely, cukpz = 
"kpzjoc = 1/2, (3 K pz = 1/3, z KPZ = 3/2. It is also 
remarkable that, in contrast to the no-overhang and small 
slopes case, Aq plays no role in the determination of the 
universality class. 

The system size of the simulations span a wide range 
of length scales, from L = 0.5 to L = 500, always en- 
suring that the continuum limit has been achieved. For 
the largest sizes, where KPZ scaling is found, we have 
used At = 5 • 10 -3 , and spatial cutoffs l m in = 5 • 10 _5 i 
and Imax = 5 • 10~ 3 L. With these values, the number 
of points composing our simulated interfaces range typ- 
ically between 100 < N < 1000. Simulation times have 
reached t = 7000 in order to ensure saturation to a steady 
state. Equation parameters have been taken of order 
unity: A 6 [0,1], A 1 £ [0,1] and A n G [0,3]. Ensem- 
ble convergence is ensured by taking, in most cases, 200 
samples or more. 

Roughness measurements are done using our diagonal- 
fitting algorithm, although checks have been performed 
with the usual Monge-gauge definition, with no change 
in the exponents. We should remark that the simulation 
algorithm is extremely robust and does not suffer from 
instabilities. 

Fig. Q] shows three typical interface profiles for L = 
75, for the same times (t = 3000), once saturation is 
achieved, for different values of the parameters. Note 
the relatively large slope values and the presence of over- 
hangs. As intuitively expected, comparison between the 
orange and the blue curves illustrates the increased abun- 
dance of overhangs for larger noise intensities. However, 
changing from the latter to the red interface, we are led 
to expecting little influence of the value of Aq on the 
morphological properties. In turn, the green line corre- 
sponds to a parameter condition that differs from that of 
the blue line in the value of A\. As seen in Fig. [2l this 
also seems to have little effect in the scaling behavior. 

The time evolution of the global interface width is 
shown in the left panel of Fig. [2] from t = 0.1 to t = 3000 
for L = 75, with 200 samples and A n = 3. For all shown 
cases, prior to saturation there is a scaling regime char- 
acterized by a exponent value /? — 0.33 ± 0.01 as given 
by a fit for the time range indicated in the figure. The 
overlap of the three curves is remarkable, this scaling be- 
havior appearing to be quite independent on parameter 
Ao being non-zero. Notice that equation ([8]) corresponds 
to non-equilibrium (irreversible) growth even if Aq = 0, 
due to removal of self-intersections. The behavior shown 
in Fig. [2] provides evidence on our claim that Aq is irrel- 
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FIG. 1. Typical interface profiles for L — 75 at saturation 
(t = 3000) for different values of the parameters as indicated. 
The curves have been artificially offset in the vertical direction 
in order to allow comparison. All units are arbitrary. 



evant to the scaling behavior. We reach a similar con- 
clusion from the measurements of additional exponents, 
although the results shown below will mostly consider 
Ao ^ conditions. 

Additionally, the value of A\ is seen in Fig. [2] to have 
little effect in the strong coupling KPZ scaling behav- 
ior. However, for smaller system sizes or smaller noise 
amplitudes, the observed value of the growth exponent 
(3 ~ 0.25 is closer to that of the Edwards- Wilkinson uni- 
versality class, see the right panel in Fig. [21 suggesting the 
system is in a prcasymptotic transient for such param- 
eter conditions, as expected from results in the Mongc 
gauge [3fij . compare both panels of our Fig. [2] with e.g. 
Figs. 2 and 4 in [H, H3|, where accurate numerical simu- 
lations of the standard KPZ equation are reported using 
a pseudospectral method. 

Regarding the value of the roughness exponent, in- 
deed Fig. [3] also shows the approach to an asymptotic 
value compatible with the KPZ universality class for suf- 
ficiently large system sizes L. Thus, the global width is 
plotted in this figure as a function of L for fixed A = 1 , 
Ai = 0.1, and A n = 3. In the figure, we estimate 
a = 0.51 ± 0.03 from a fit of the behavior for the largest 
system sizes 75 < L < 200. In general, the transition 
from non-KPZ prcasymptotic (small L, A n values, weak 
coupling regime) to asymptotic KPZ (large L, A n val- 
ues, strong coupling regime) scaling behavior is not par- 
ticularly correlated with an increasing dynamical role of 
self-intersection removal. The simulations show that the 
length of interface that is removed per unit time reaches 
a constant non-negligible value for generic parameter and 
system size conditions. 

The large slopes and overhangs seen in Fig. [T] suggest 
a study of the statistics of the local angle 9 shown by the 
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FIG. 2. Time evolution of the global interface width for L = 
75 and different parameter choices. Left panel: for A n = 
3, the scaling behavior does not depend appreciably on the 
values of Ao and Ai. The straight line is a fit to the data 
within the corresponding time interval, its slope being close 
to the KPZ behavior /?kpz = 1/3. Right panel: comparison 
between A„ = 3 case shown in the left panel and A n = 1. 
The straight lines are obtained from fits to the data within 
the corresponding time intervals, having slopes that are close 
to the KPZ and EW behaviors, /3kpz = 1/3 and /few = 1/4, 
respectively. All units are arbitrary. 



interface with respect to the global horizontal direction. 
The left panel of Fig. [4] shows the angle histogram (i.e., 
angle probability distribution) for L — 75 and different 
values of the parameters, once saturation was achieved. 
A fit is shown to an exponential distribution of the form 
C\ exp(— C2 1 6* I x ) • Again, A does not seem to play any 
role in this behavior, and the three curves overlap. The 
right panel shows the evolution of this histogram as L 
increases. Fits to the same exponential form are also dis- 
played. The corresponding values of \ increase with L, 
stabilizing to a value close to 2, i.e., the distribution be- 
comes Gaussian. In any case, the symmetric exponential 
distribution (as opposed to a power-law one) centered 
around the origin implies that the system has a charac- 
teristic slope that is horizontal. Of course, the probability 
distribution can not be truly Gaussian, since the angle 
domain is bounded in [— 7r, vr]. 
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FIG. 3. Global interface width as a function of system size. 
For comparison, the straight line is obtained from a fit of the 
behavior for 75 < L < 200 and has a slope corresponding to 
a = 0.51, compatible with q?kpz = 1/2. In all cases Aq = 1, 
Ai = 0.1, and A n = 3. All units are arbitrary. 



The change in the local morphological properties for in- 
creasing size L goes beyond the angle histogram. Thus, 
the local roughness exponent ai oc is close to 1 /2 for small 
noise amplitudes A n < 1. If the noise amplitude is large 
{A n = 3), for small systems (L < 50) this exponent takes 
an effective value that can raise even above 0.6, see left 
panel on Fig. [5] where the local interface width at satu- 
ration, w(l), is shown as a function of the local scale, for 
different system sizes L and fixed Ao = 1, A± = 0.1, and 
A n = 3. Moreover, the non-complete overlap of e.g. the 
curves for L = 25, 50, is reminiscent of anomalous scal- 
ing behavior. However, for larger values of L > 75 the 
local roughness exponent converges to an L-independent 
value a^c =0.51 ±0.01, namely, close within error bars 
to the 1/2 value it takes both for the standard EW and 
KPZ equations in ID. Moreover, for these system sizes 
curves overlap, a fact that indicates that asymptotic scal- 
ing is not anomalous. The right panel in the same fig- 
ure shows how the local width w(l, t) curves evolves with 
time for our largest system size, L = 200. The curves 
give the behavior of the local width with I for different 
times. Apart from times that are comparable with the 
discretization unit, for which no proper correlations have 
built up, curves for increasing times overlap in the form 
that again is expected for FV scaling. As time evolves, 
the lateral correlation length £(t) increases, finally reach- 
ing a steady state value at which a\ oc ~ 0.5 is obtained 
from a fit of the large I behavior, compatible with the 
value of the global roughness exponent a found above. 

Actually, we can employ the local morphological be- 
havior just considered in order to measure the dynamic 
exponent z, as an alternative check to assess the asymp- 
totic behavior of our system. As mentioned above, within 
the Family- Vicsek Ansatz, the saturation time t s ~ L z , a 
relation which is also fulfilled for local scales, t s {l) ~ l z . 
This implies that, if time is measured in units of V and 
the local interface width is measured in units of l a , all 
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FIG. 4. Left panel: angle histograms are shown at saturation 
for different choices of Ao and Ai, L = 75, and A„ — 3, 
obtaining very similar results. The solid line is a fit to a shape 
Ci exp(-C 2 |6>| x ) with x = 1-54. Right panel: for fixed A = 
1, A\ = 0.1 and increasing size, the angle histograms can be 
seen to evolve towards a Gaussian profile. Solid lines are fits 
similar to that in the left panel with % = 1-44, 1.54, 1.84, 1.99 
for increasing L. 



the w(l,t) curves should collapse. This is shown in Fig. 
[6] for the L = 200 curves shown on the right panel of Fig. 
E As usual (H, the values of a and z are chosen so that 
the overlap between the curves is maximum. This occurs 
for a = 0.52 ± 0.2 and z = 1.4 ± 0.1, again in agreement 
with the KPZ values for the roughness and the dynamic 
exponents, respectively. 



VI. CONCLUSIONS AND OUTLOOK 

We have characterized the scaling properties of a con- 
tinuum stochastic model of one-dimensional interfaces, 
given by Eq. (|9|). The model captures the main physical 
mechanisms believed to characterize non-conserved, irre- 
versible surface growth, as described in particular within 
the KPZ approach to non-conserved interface dynamics, 
namely, growth along the local normal direction with a 
rate that is affected by local curvature and by fluctua- 
tions in growth events, but without assumptions on small 
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FIG. 5. Left panel: local interface width as a function of scale, 
w(l), at saturation for different system sizes as indicated. In 
all cases, Ao — 1, A\ = 0.1, and A n = 3. The two groups 
of curves (L < 75 and L > 75) have been artificially offset 
for clarity. Notice the overlap between the different curves 
for L > 100, implying that scaling is FV. The small L < 
100 behavior is slightly anomalous, and characterized by a 
different value of a\ oc - The solid lines are obtained from fits 
for the smallest and largest values of L in the corresponding 
regions of I. Their slopes provide the value ai oc in each case. 
Right panel: local interface width w(l,t) vs scale for L — 200 
and logarithmically equispaced times between t = 0.1 and 
t = 200, bottom to top. The straight line with slope ai oc = 0.5 
is obtained from a fit of the behavior for the top curve within 
the corresponding range of /. All units are arbitrary. 



slopes or absence of overhangs. Our numerical simula- 
tions are tailored to avoid these simplifications, and lead 
to a number of results. In spite of the possibility for the 
surface to develop arbitrarily large slopes, a fact that is 
frequently associated with anomalous scaling, the kinetic 
roughening properties fulfill, rather, the Family- Vicsek 
Ansatz. There are predictions on the non-occurrence of 
non-trivial anomalous scaling at the steady state of lo- 
cal growth equations [4l[. In principle, these predictions 
are done within a small slope approximation so that our 
present result supports their validity beyond such con- 
straint for the specific example that we have studied. 
This suggests in turn the interest of assessing the cor- 
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FIG. 6. Data collapse of the local interface width w(l, t) 
curves of the right panel on Fig. [5] [L = 200). Collapse is 
achieved for a = 0.52 and z — 1.4. All units are arbitrary. 

responding situation for other continuum models beyond 
the small slope approximation. 

Regarding the specific type of scaling that we find, 
for system sizes and/or noise amplitudes that are small, 
the system is in a prcasymptotic state within which the 
growth exponent is close to that of the EW universality 
class, while the roughness exponent is close to but not 
necessarily equals the corresponding EW value. In any 
case, for sufficiently large system sizes and/or noise am- 
plitudes, Eq. ([9]) turns out to be in the KPZ universality 
class. Naturally, this result agrees with general consider- 
ations based on symmetries and conservation laws [l|, @] , 
but we find it to be non-trivial. To the best of our knowl- 
edge, it is found for the first time for a SPDE that can 
access height configurations that are not constrained by 
the small slope and no-overhang approximations. Note 
moreover that, as a result of the very formulation of Eq. 
©, the role of noise in this system is very far from the 
standard KPZ case. Thus, in the latter noise is decorre- 
lated only in the x-axis. Adoption of that type of noise in 
our system would amount to assuming that all interface 
points which share the same x-coordinate are correlated, 
disregarding their distance in the vertical direction. In 
our case, noise is decorrelated both in the horizontal and 
vertical directions, as seen in Eq. (0. Thus, its naive spa- 
tial scaling dimension should be 2 rather than 1. Even a 
power counting argument [l[ on the KPZ equation with 
this modification leads to a different result, since it pro- 
vides consistent scaling exponents which differ from the 
KPZ universality class. Thus, our present numerical re- 
sults are not evident from this point of view. 

In our simulations, we have seen moreover that our 
extra condition on self-intersection removal plays a non- 
negligible role, in particular driving the system out of 
equilibrium even when Ao = 0. Note that the equilibrium 
problem associated with Eq. ([9]) is well understood (42| . 
In our case, intersection removal is the price we have to 
pay in order to have mathematically a simply-connected 
interface. Physically, it might be thought of as playing 
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a similar simplifying role as the condition of a phantom 
membrane plays with respect to more realistic polymer- 
ized membranes [13] • Still, we believe the ensuing model 
features non-trivial physics. Perhaps this can be appre- 
ciated better if we think of still more realistic models of 
growth, of which Eq. (O and similar equations are sim- 
plified versions. In those, the dynamics of the interface 
is coupled to that of a physical field, such as e.g. the 
concentration of aggregating units in the dilute phase, in 
what mathematically is described as a (stochastic) mov- 
ing boundary problem, see e.g. [U |43|. For parame- 
ter conditions in which kinetic roughening occurs, voids 
and bubbles can indeed be createdfleading to a multi- 
ply connected interface), see e.g. [30|, but their role for 
the morphological evolution seems to be quite marginal 



31|. Incidentally, self-intersection removal seems to be 
a relevant stabilizing mechanism in face of morphologi- 
cal instabilities, as in ox-bow formation in meandering 
rivers [3, EH . In stable-growth moving boundary prob- 
lems, the main role in the evolution is played, rather, 
by the envelope of the so-called active growth zone, that 
corresponds in our model to the simply connected inter- 
face that we keep track of after self-intersection removal. 
From this point of view, our choice does not seem too 
distant from the standard procedure by which the full 
interface dynamics of e.g. discrete growth models that 
lead to bubbles and overhangs due to e.g. bulk vacancies 
[46| is traded for that of its single valued approximation. 

An interesting future extension of the present work is 
the study of the full Eq. (JSj> , namely, for a non-zero 



coefficient, that would allow us to study surface diffu- 
sion within the present framework, with potential im- 
plications in the context of thin film growth by MBE 
[10( . Given the slow dynamics associated with relaxation 
by surface diffusion, this would require improvements on 
the more computationally expensive ingredients of our 
adaptive numerical scheme, conspicuously including the 
self-intersection removal algorithm. 

Moreover, given the geometrically intrinsic nature of 
the dynamics, the numerical approach, and the mea- 
surement techniques developed in the present work — by 
which nonlinearities and noise correlations, which would 
appear in the Monge gauge, come out naturally — , ad- 
ditional interesting applications would be the study of 
growth in a circular geometry (see a recent discussion 
in p7j and references therein) and the extension to two- 
dimensional interfaces. In both cases there are still in- 
teresting open theoretical and experimental challenges. 
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